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Method for evaluating the operating conditions of a machine or an 
installation 

The invention relates to a method for evaluating the operating 
5 conditions of a machine or an installation. 

Until now the operating conditions of a machine or an installation 
were evaluated using visual extrapolations of a, for example, 
critical measurement value pattern and/or subjective evaluations of 
10 the impact on parameters of the machine or the installation, in 
order to estimate the operating conditions of a machine or an 
installation and to respond accordingly by changing or otherwise 
influencing the parameter according to the evaluation, for example 
by predetermining a target value for parameters . 

15 

In this process adjustment /regress ion functions of a database were 
adjusted and an optimization process was carried out by iterative 
selection of curve functions, with the maximum correlation 
coefficient. The curvature pattern of such an adjusted curve does 

20 not necessarily have to correspond to that of the database. The 
correlation coefficient r (maximum value = 1, minimum value = -1) 
can only be used as an adjustment quality criterion under certain 
conditions, as this value depends not only on the adjustment quality 
of the curve function used but also on the gradient of the curve 

25 function used. If the gradient, for example, of a linear adjustment 
tends toward zero, r also follows this trend, regardless of the 
scatter of the individual curve points. 

This means that r cannot be used as a measure of quality for an 
30 extrapolation. 

Generation tools for evaluating the operating conditions of a 
machine or an installation have to satisfy certain minimum 
requirements, so that an evaluation of the extrapolation result is 
35 possible in respect of 
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• Predictive reliability 

• Variables influencing predictive dependability 

• Traceability of the prediction. 

5 Naturally any prediction is subject to uncertainty and its measure 
of quality for assessing forecast alarms /exceeded limit values and 
as a basis for decision for resulting, e.g. automated, actions is 
extremely important. In addition, in the case of a cyclically 
generated method for evaluating the operating conditions of a 

10 machine or an installation in a specific, e.g. measurement context, 
the change in a measure of quality over time can also be seen as a 
trend and can therefore also offer additional conclusions about the 
time of occurrence of the forecast event, e.g. the exceeding of a 
limit value, in order to ensure the operational dependability of the 

15 machine or the installation. 

Estimating the reliability of an extrapolation is of central 
significance, as evaluations of damage symptoms, the exceeding of 
limit values and operational optimization, e.g. replacement of parts 
20 is not really possible without knowledge of the predictive 
reliability, in other words the measure of quality. 

An evaluation of the operating conditions of a machine or an 
installation, in which a trend in the further changes of parameters 

25 is analyzed, is very important in respect of problem solving in 
tasks such as operation monitoring or system analysis. 
These methods are used wherever extrapolations into unknown value 
ranges from known laws are of significance, in other words 
essentially when predicting events (identification of future 

30 alarms /exceeded load values) and damage prevention (identification 
of future damage) . 
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As with measurement values, these are practically worthless without 
details of measurement tolerance, a variable must be associated with 
an extrapolation, allowing conclusions to be drawn about predictive 
quality (measure of quality) . Without this information extrapolation 
5 results cannot be properly evaluated and decisions on follow-up 

actions (e.g. process interventions) are extremely unreliable and in 
some circumstances counterproductive. 

The object of the invention is therefore to eliminate this problem. 

10 

This object is achieved by means of a method according to Claim 1. 

Further advantageous embodiments are listed in the subclaims. 

15 Optimal adjustment of the extrapolation function in respect of the 
database is achieved not by a maximum correlation coefficient r but 
using a measure of quality K, which is independent of the gradient 
of a curve of adjustment in respect of a database. 

20 The measure of quality K is calculated from at least two different 
subcriteria and this provides conclusions about the predictive 
quality of the extrapolation. The existence of K allows at least two 
extrapolation modes : 

25 • Predetermination of a constant measure of quality K 

In this case the extent of the extrapolation, i.e. the distance 
between the last value in the database and an end point of the 
extrapolation, is variable. 



30 
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This procedure is important when suppressing incorrect or 
insignificant alarms. 

Here an operator or technical installation can be informed of the 
period for which no alarm will occur, i.e. a specific operational 
5 dependability is specified. 

• Predetermination of the extent of extrapolation (e.g. when a 
limit value is exceeded) 

10 In this case K is variable. The change to K over time can in turn be 
used as a basis for a trend analysis, in order to estimate the 
probable time of occurrence of the event when K = 1. 

Trend analyses to evaluate the operating conditions of a machine or 
15 an installation can generally extrapolate the change in two 
correlated variables into unknown value ranges . 
They can be applied both to changes over time 

y± = f (tj ; i = l...n 

Example P e i = f(t); P e i = generator output 

20 and formally time- independent associations 

yii = f(y2i); i = l...n, yi synchronous with y 2 

Example P e i = f (fuel mass flow) . 

The use of trend analyses however implicitly assumes that conditions 
25 of laws that can be derived from the past for changes in two 

correlated variables also continue in the future and unknown value 
associations can be generated in this way. For meaningful tend 
analyses, this condition requires that these can only be generated 
with patterns which were generated under identical (standardized) 
30 basic conditions - physical environmental 
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conditions, fuel, output, burner operating mode, etc. Value 
associations standardized in this way can be supplied either by data 
selection using logical conditions (discrete standardization) or 
using a process simulation model (analytical standardization) , which 
5 can for example calculate all environment-dependent values back to 
ISO conditions. 



The figures show 



10 Figure 1 a database and extrapolation curve 

Figure 2 curves to determine AI 

Figure 3 segmentation of a database 

Figure 4 a pattern over time of K 

Figure 5 a turbine with components 

15 Figure 6 a simplified diagram of a longitudinal section of a 

turbine and 

Figure 7 a turbine blade as a component . 



Basic variables for calculating measure of quality K 

20 

The measure of quality K is a measure of the quality of an 
extrapolation. It links a number of error influences and is 
therefore a function of at least two variables, e.g. V, AI. 
In this case these are for example a ratio V, an x uncertainty AI, a 
25 continuity S and a time constant C: 



K = f (V, AI, S, C) . 



These variables are for example linked so that the value repertoire 
30 of K is within a standardized range, e.g. 0% to 100%. 
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• Ratio V of the data range to the extrapolation range 

There are n measurement points (xi, yi) to (x n , y n ) available for the 
extrapolation. However not all the available measurement points 15 
5 of a data repertoire 54 (Fig. 3) have to be used. 

The ratio V for example inserts the x value range {xx-xj of the 
database 6 used for the extrapolation into the ratio of the distance 
Xx to the x position of the extrapolation point x s used (Fig. 1), 
i.e. the extrapolation range 9, into a ratio, which for example aims 
10 as a maximum toward the value 1, i.e. the database 6 is equal to the 
extrapolation range and as a minimum toward the value 0, i.e. the 
database 6 is very small compared with the extrapolation interval. 
X s is for example a value, at which a limit value (Y s = limit value) 
is exceeded. 

15 Figure 1 also shows an adjustment function 12 in respect of the 
database 6, which also gives the extrapolation to the x point x s . 

V is for example calculated as follows: 

20 V = (x n - x x ) / (x s - x x ) . 

If V essentially contributes to a reduction in the measure of 
quality K f the database of the extrapolation can gradually be 
extended. 

25 

• Uncertainty AI of the adjustment curve in the x direction 

Every curve adjustment (condition for extrapolation) naturally also 
has an uncertainty (confidence range) of the curve parameter 
30 calculated. 

In order to determine the uncertainty of the gradient, in Figure 2 
for example a linear regression function 21 is used. 
Generally all functions, which can be transferred to linear 
structures, can be used as extrapolation functions, i.e. 



35 
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Linear function (master function) -> y = a 0 + a x *x 

Potency function -> In y = lna 0 + a^lnx 

Logarithmic function -> y = a 0 + ai *ln x 

Exponential function -> In y - In a 0 + a x *x 
etc . 



The uncertainty AI is then calculated as follows. 

Figure 2 shows an example of a database 6 with an extrapolation 

curve 21 and with further curves 24, 27 for calculating AI . 

The uncertainty AI can be illustrated by rotating the adjustment 
curve 21 (shown in a linear manner) about a point of rotation (x, y) . 
The mean values x and y are for example determined by geometric or 
arithmetic means. The point (x,y) does not necessarily have to lie on 
the linear adjustment curve 21 (extrapolation curve) . 
The uncertainty of a gradient of the linear curve 21 corresponds to 
an angular rotation, with which the gradient m changes by ±Am and 
therefore an x uncertainty AI of the adjustment curve, in the x 
direction for example (Fig. 2). 

Varying the gradient results in two further curves 24, 27, each of 
which has a point of intersection 30, 33 with a parallel 18 (ys = 
constant), which corresponds to a limit value. For each point of 
intersection 30, 33, there is a corresponding x value 1^ and I min , 
whereby > and AI = 1^ - i mln . 

This x uncertainty correlates to an intersection angle a, which is 
determined from the angle between the parallel 18 and the curve at 
the point of intersection of the extrapolation curve 21 and the 
parallel 18. The x uncertainty AI increases as the intersection 
angle a reduces. 

If there is no limit value or no limit value is exceeded with an 
extrapolation, the above applies to the predetermined x end point 
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(x s ) of the extrapolation and an imaginary parallel to the x axis, 
which runs through the point {x s/ y s ) . 

As the extrapolation function always exists analytically, the 
intersection angle a between a horizontal or alarm function and the 
extrapolation function can be calculated using the first derivative. 

The approach to calculating the gradient uncertainty ±Am is for 
example as follows: 

A gradient b is determined by means of a confidence interval. Such a 
method is known from Kreyszig, Erwin: "Statistische Methoden und 
ihre Anwendungen" [Statistical methods and their applications] , 
published by Vandenhoeck und Ruprecht, Gottingen, page 270. 

The following procedure is used to determine the gradient b: 

a) A database 6 is established. The database 6 comprises n 
correlated x and y values (Fig. 2) . 

b) Calculation of x and y of the database 6 and the variable 



c) 



Calculation 




(i = 1 ... n) 



d) 



Calculation 




e) 



Calculation 




(i = 1 ... 



n) 
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g) Calculation of a = (n -l)(y* -f/g 2 ) 

The following equation is obtained for the regression line 21 

y = y + b(x-x). 

Next a confidence factor of for example 95% is predetermined, 
i.e. y = 0.95, from which a variable F (c) is calculated: 

h) F(c) = 1(1 + y) = 0.975 

i) With F(c) = 0.975 and n-2 (n = number of measurement values) 
degrees of freedom, the t-distribution (student distribution) 
gives a value c 

(0.975 corresponds to the integral of the t-distribution to 
the point x = c) . 

Am is obtained from 
j) C ^~ a 

This gives an uncertainty of the gradient m: 
b-Am<m<b + Am. 

This calculates the x uncertainty Al of the gradient m in the 
point of rotation (Fig. 2) depending on the scatter variables of 
the current database 6. 

As a result Al = f (Am, a) can be calculated using radiation or 
goniometry . 

Al can be standardized similarly. 
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• Continuity S of the y values in the database 

Segmentation of a database 6 provides prior conclusions about the 
current direction of curvature of the database 6 used. 
5 The database 6 is a true or false subset of a data repertoire 54 

comprising all available measurement points 15. 

Segmentation means that measurement points, i.e. numeric tuples 
(x,y) of the database 6 are divided into at least three segments 
45, 48, 52. A linear adjustment curve 36, 39, 42 is determined for 
10 each segment 44, 48, 52. 

The direction of curvature of the database 6 is determined by 
creating segment means gl, g2 and g3 and calculating the second 
numerical derivative. 

15 With this method it can be assumed that the adjustment curve 

(extrapolation curve) is calculated cyclically from a floating 
database 6, i.e. a data window 3 (Fig. 3) of constant or variable 
length is displaced after each completed extrapolation cycle in 
the direction of increasing variables by a predetermined interval 

20 Ax (independent variable) in each instance. 

In practice the displacement Ax of the data window 3, i.e. the 
database 6, can be effected to the maximum in the time cycle, with 
which new measurement values 15 are generated (e.g. Ax = 5 s) . 
The concept of continuity S in this context is not to be 

25 understood as a mathematical definition but as a measure of the 

change in the point pattern of the database 6 in relation to the 
last completed extrapolation stage. 

In order to be able to make a statement about the continuity of 
the Y value change and therefore implicitly a statement about the 

30 extrapolability of the change, the data window 3 (Fig. 3) is 

divided into at least three segments 45, 48, 52. In an 
intermediate stage, a mean and a linear adjustment function y l7 
Y2i Y3 with the gradients Ci, c 2 , c 3 are calculated for each 
segment 45, 48, 52 (Fig. 3). 

35 If the means of the three regression lines 36, 39, 42 are 
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designated as g lf g 2 {middle segment) and g 3 , the current direction 
of curvature of the gradient pattern can be determined with the 
numerical curvature measurement p 
p = g x - 2*g 2 + g 3 . 

5 

In the example in Fig. 3 q is negative. This means a right 
curvature . 

With p = 0 the three means lie on a straight line. From a curve 
repertoire, that curve type is adjusted iteratively in respect of 
10 the value range of the entire current database. 

Adjustment functions (regression functions) of the curve repertoire 
for the extrapolation curve for data correlations have to satisfy 
the monotonous pattern condition, as non-monotonous functions can 
15 result in significant extrapolation uncertainties. Generally all 

functions which can be transferred to linear structures of the form 
y = a 0 + ai*x can be used, i.e. 



Linear function (master function) 
20 Potency function 

Logar i thmi c f unc t ion 
Exponent i a 1 f unc t i on 
etc. 



-> f (x) = y = a 0 + ai*x 

-> f (x) = In y = lna 0 + a^lnx 

-> f (x) = y = a 0 + ai*ln x 

-> f (x) = In y = In a 0 + ai*x 



25 The curve type selected from the curve repertoire must satisfy the 
following conditions: 

the direction of curvature of the extrapolation curve (regression 
lines) must correspond to that of p; 

30 the quotient Q k of numerator = mean (weighted if necessary) of the 
distance squares between measurement values and extrapolation curve, 
and denominator = mean square of the y value range of the 
extrapolation curve in the area of the data window (here for example 
y 2 mitt k = [ (ymaXk + ymink)/2] 2 , where ymax is the maximum Y value and 

35 ymin correspondingly is the minimum Y value of the kth curve) and is 
used for standardization) must be minimal; 
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Q k =/(*)= 



where 

k is a numerator of the available extrapolation curve types (curve 
5 repertoire) 

Yi(xi) is the measurement value at point x± 

f k (Xi) is the function value of the kth extrapolation curve type at 
point Xi 

Wi is a weighting factor for each individual measurement value or for 
10 all measurement values of a segment. 

The continuity S is evaluated as follows: 

Comparison of the three straight line gradients d, C 2 , C 3 with the 
gradients of the extrapolation curve at the respective middle 

15 positions (X s i, X S2 , X S3 ) of the three segments of the data window. A 
different weighting of the three gradient differences is possible 
and expedient here. In this way the more current values in the last 
segment can be evaluated more stringently, in order to identify a 
change in the curve pattern more rapidly. 

20 The three gradient differences are a measure of the continuity of 
change in the curve pattern. 

S is calculated as follows: 

Oi to 0 3 are the three gradients of the selected kth adjustment curve 
25 36, 39, 42 in relation to the half segment width in each instance, 
Ci to C 3 the gradients of the linear segment adjustments (Fig. 3). 

The following therefore applies: 



30 




(segment areas) 



with yi: weighting factors l...n. 
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The value range of S is as follows: 0 (abs. continuity), i.e. the 
gradients of the extrapolation function in the segment centers are 
identical to the regression lines of the individual segments 45, 48, 
52 to + oo. i.e. no gradient correspondence. 

5 

• Time constancy C of the extrapolations 

The time variance of a sequence of extrapolations is an additional 
and important indicator of predictive dependability, as high 
10 variance totals for example of the measure of quality can testify to 
non-stable extrapolation conditions and therefore questionable 
significance. 

The time constancy of an extrapolation can be requested when a fixed 
measure of quality is predetermined with a variable extent of 
15 extrapolation (X s ) or when the extent of extrapolation is fixed - 
e.g. X s can correspond to an exceeded fixed limit value - via the 
then variable confidence factor. 

Figure 4 shows an example of the time pattern of the measure of 
quality K(t) with different values 58. For this value pattern an 
20 adjustment function q is determined from a curve repertoire. 

The respective variables are calculated iteratively, e.g. using the 
Newton method. 

Variances are monitored over time using a regression fit (e.g. by a 
polynomial 1 or higher order, with the order dependent on n) in 

25 respect of the last n (e.g. n = 10) time values for extent of 
extrapolation or confidence factor. In this way values for the 
variables to be tested, which change in a linear manner over time, 
can be evaluated with precision in their variability. 
The evaluation of the time variance of the extent of extrapolation 

30 is calculated as follows: 
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i = number of iterations. 

5 Here too weighting factors Yi can evaluate the most current values as 
different from/higher than the oldest. 
Here q 2 mitt k corresponds to y 2 mitt k . 

The value range for C is between C=0 (no temporal scatter of 
variables, all variable points are on the regression line) and + «> 
10 (variable sequence degraded) . 

• Linking of V, Al norm , S and C to measure of quality K 

The measure of quality K links the variables V, Al norTn/ S and C for 
15 example as follows: 

„ V*A7 

A — . 

S*C 

The value range of K is between 0 (extrapolation distorted) and + <~ 
20 (extrapolation exact) . 

K can for example be standardized by means of a non-linear equation: 

Knorm ^ & 

25 The installation is for example a rotor of a gas turbine. By 

determining K of the frequency of rotor parameter, it is ascertained 
that in a specific period (xs-xn) a limit value will be exceeded 
without permission. The frequency is therefore controlled downwards. 
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Figure 5 shows a diagram of a longitudinal section of a gas turbine 
1. 

A gas turbine 1 is selected as an example of a component for a 
5 machine, in which parameters such as temperature, vibration, 
electrical power or other parameters are measured. 
A compressor 7, combustion chamber 10 and turbine component 13 are 
arranged one behind the other along a shaft 4. The turbine component 
13 has a hot gas duct 16. Gas turbine blades 20 are arranged in the 
10 hot gas duct 16. Guide blade and runner blade limits are arranged 
alternately one after the other. 

The gas turbine blades 20 are for example cooled by means of a 
combined air and/or steam cooling system. For this compressor air is 
for example taken from the compressor 7 and fed to the gas turbine 
15 blades 22 via an air supply system 23. Steam is also fed to the gas 
turbine blades 20 for example via a steam supply system 26. 

Figure 6 shows a simplified diagram of the longitudinal section of a 
turbine component 13. The turbine component 13 has a shaft 4, which 
20 extends along an axis of rotation 41. 

The turbine component 13 also has an inflow area 49, a blade area 51 
and an outflow area 53 one after the other along the axis of 
rotation 41. 

Rotatable runner blades 20 and fixed guide blades 20 are arranged in 
25 the blade area 51. The runner blades 20 are attached to the shaft 4, 
while the guide blades 20 are arranged on a guide blade support 47 
around the shaft 4. 
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A ring-shaped flow duct for a flow medium A, e.g. hot steam, is 
formed by the shaft 4, the blade area 51 and the guide blade support 
47. The inflow area 49 used to feed in the flow medium A is limited 
in a radial direction by an inflow housing 55 arranged upstream from 
5 the guide blade support 47. 

An outflow housing 57 is arranged downstream on the guide blade 
support 47 and limits the outflow area 53 in a radial direction, 
i.e. perpendicular to the axis of rotation 41. During the operation 
of the gas turbine 1, the flow medium A flows from the inflow area 
10 49 into the blade area 51, where the flow medium works subject to 
expansion, and then leaves the gas turbine 1 via the outflow area 
53. The flow medium A is then collected in a condenser (not shown in 
more detail in Figure 5) for a steam turbine beyond the outflow 
housing 57. 

15 

As it flows through the blade area 51 the flow medium A expands and 
works on the runner blades 20, causing these to rotate. 

Figure 7 shows a perspective view of a runner blade 20, which 
20 extends along a radial axis 60. The runner blade has an attachment 
area 63, an adjacent blade platform 66 and a blade paddle area 69 
one after the other along the radial axis 60. 

A blade foot 72 is formed in the attachment area 63 and is used to 
attach the blade 20 to the shaft 4 of a gas turbine 1. The blade 
25 foot 72 is for example configured as a hammer head. 

With conventional runner blades 20 solid metal materials are used in 
all areas 63, 66, 69. The runner blades 20 can be manufactured 
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using a casting method, a forging method, a milling method or 
combinations of these. The component frequently exhibits defects 
immediately after manufacture. 



